Electron- hydrogen scattering in Faddeev-Merkuriev integral equation approach 
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Introduction 

The scattering of electrons on hydrogen atom is a fun- 
damental three-body problem in atomic physics. The 
long-range Coulomb interaction presents the major diffi- 
culty. On the other hand, it is a special kind of Coulomb 
three-body problem as it contains two identical parti- 
cles. While many studies have been carried out aiming 
at solving the Schrodinger equation using perturbative, 
close-coupling, variational or direct numerical methods 
approaches along to the Faddeev equations are relatively 
scarce. Here, by solving Faddeev-type integral equations, 
we present a general numerical method suitable for the 
treatment of elastic and inelastic processes in three-body 
Coulombic systems with two identical particles and apply 
the formalism for the electron-hydrogen system. 

For quantum mechanical three-body systems the Fad- 
deev integral equations are the fundamental equations. 
They possess connected kernels and therefore they are 
Fredholm-type integral equations of second kind. The 
Faddeev equations were derived for short range inter- 
actions and if we simply plug-in a Coulomb-like poten- 
tial they become singular. The necessary modification 
were proposed by Merkuriev Q . In Merkuriev's approach 
the Coulomb interactions were split into short-range and 
long-range parts. The long-range parts were included 
into the ,,free" Green's operators and the Faddeev pro- 
cedure were performed only with the short-range poten- 
tials. The corresponding modified Faddeev, or Faddeev- 
Merkuriev equations are mathematically well-behaved. 
They possess compact kernels even in the case of at- 
tractive Coulombic interactions. This means that the 
Faddeev-Merkuriev equations possess all the nice prop- 
erties of the original Faddeev equations. 

However, the associated three-body Coulomb Green's 
operator is not known explicitly. To circumvent the prob- 
lem the integral equations were cast into differential form 
and the appropriate boundary conditions were derived 
from the asymptotic analysis of the three-body Coulomb 
Green's operator. These modified Faddeev differential 
equations were successfully solved for various atomic 
three-body problems, including electron- hydrogen scat- 
tering up to the H(n = 3) threshold 

A characteristic property of the atomic three body 
systems is that, due to attractive Coulomb interactions, 



they have infinitely many two-body channels. If the to- 
tal energy of the system increases more and more chan- 
nels open up. The differential equation approach needs 
boundary conditions for each channels, and becomes in- 
tractable if the energy increases beyond a limit. Integral 
equations do not need boundary conditions, this informa- 
tion is incorporated in the Green's operators. They need 
initial conditions, which are much simpler. Therefore an 
integral equation approach to the three-body Coulomb 
problem would be very useful, it could provide an unified 
description of the scattering and reactions processes for 
all energies. 

In the past few years we have developed a new ap- 
proach to the three-body Coulomb problem. Faddeev- 
type integral equations were solved by using the 
Coulomb-Sturmian separable expansion method. The 
approach was developed first for solving the nuclear 
three-body scattering problem with repulsive Coulomb 
interactions ||, which has been adapted recently for 
atomic systems with attractive Coulomb interactions 
M. The basic concept in this method is a ,, three- 
potential" picture, where the S matrix is given in 
three terms. In this approach we solve the Faddeev- 
Merkuriev integral equations such that the associated 
three-body Coulomb Green's operator is calculated by 
an independent Lippmann-Schwinger-type integral equa- 
tion. This Lippmann-Schwinger integral equation con- 
tains the channel-distorted Coulomb Green's operator, 
which can be calculated as a contour integral of two-body 
Coulomb Green's operators. The method were tested 
in positron-hydrogen scattering for energies up to the 
H(n = 2) — Ps(n = 2) gap Q, and good agreements 
with the configuration-space solution of the Faddeev- 
Merkuriev equations were found. 

In this paper we apply this formalism for the electron- 
hydrogen scattering problem. In Sec. I we briefly describe 
the Faddeev-Merkuriev integral equations, the details are 
given in Ref. [Q. However, the fact that in the electron- 
hydrogen system we have to deal with identical particles 
requires some additional considerations: the symmetry 
simplifies the numerical procedure. In Sec. II the inte- 
gral equations are solved by the Coulomb-Sturmian sep- 
arable expansion method. In Sec. Ill we show some test 
calculations up to the H(n — 4) threshold with total an- 
gular momenta L = and L = 1. Finally, we draw some 
conclusions. 
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I. FADDEEV-MERKURIEV INTEGRAL 
EQUATIONS FOR THE e~ + H SYSTEM 

In the e~ + H system the two electrons are identical. 
Let us denote them by 1 and 2, and the non-identical 
proton by 3. The Hamiltonian is given by 



H = H { > 



+ + v- 



3 i 



(1) 



where H° is the three-body kinetic energy operator and 
i>„ denotes the Coulomb interaction in the subsystem a. 
We use the usual configuration-space Jacobi coordinates 
x a and y a , where x a is the coordinate between the pair 
(/3, 7) and y a is the coordinate between the particle a and 
the center of mass of the pair (/?, 7). Thus the potential 
v^, the interaction of the pair (/?, 7), appears as v^(x a ). 
The Hamiltonian (0) is defined in the three-body Hilbert 
space. So, the two-body potential operators are formally 
embedded in the three-body Hilbert space, 



\x)l 



(2) 



where \ y is a unit operator in the two-body Hilbert space 
associated with the y coordinate. 

The role of a Coulomb potential in a three-body sys- 
tem is twofold. In one hand, it acts like a long-range po- 
tential since it modifies the asymptotic motion. On the 
other hand, however, it acts like a short-range potential, 
since it correlates strongly the particles and may support 
bound states. Merkuriev introduced a separation of the 
three-body configuration space into different asymptotic 
regions Q . The two-body asymptotic region is defined 
as a part of the three-body configuration space where the 
conditions 



\x\ <x a (l + \y\/y a ) 1 ^, 



(3) 



with Xo,yo > and v > 2, are satisfied. Merkuriev 
proposed to split the Coulomb interaction in the three- 
body configuration space into short-range and long-range 
terms 



V C = y (s) 



.(0 



(4) 



where the superscripts s and / indicate the short- and 
long-range attributes, respectively. The splitting is car- 
ried out with the help of a splitting function £, 



v {s) {x,y) = v c (x)({x,y), 
v®(x,y) = v c (x)[l-((x,y)} 



(5a) 
(5b) 



The function £ vanishes asymptotically within the three- 
body sector, where x ~ y — * 00, and approaches one in 
the two-body asymptotic region f2, where x « y — > 00. 
Consequently in the three-body sector vanishes and 
v^' approaches v c . In practice usually the functional 
form 

C(x, y) = 2/ {I + exp [{x/x y /{I + y/y )}} , (6) 



is used. 

In the Hamiltonian (0) the Coulomb potential v§ , the 
interaction between the two electrons, is repulsive, and 
does not support bound states. Consequently, there are 
no two-body channels associated with this fragmentation. 
Therefore the entire v§ can be considered as long-range 
potential. Then, the long-range Hamiltonian is defined 
as 



£T« = H 







,(0 



,(0 



(7) 



and the three-body Hamiltonian takes the form 

H = H® +v[ s) +v ( 2 s \ (8) 



So, the Hamiltonian (||) appears formally as a three- 
body Hamiltonian with two short-range potentials. The 
bound-state wave function |<3>) satisfies the homogeneous 
Lippmann-Schwinger integral equation 

I*) = G« U s) + 4 S) ] I*) = G®v[ a) \9) + G«i4 s) |*), 

(9) 

where G^(z) = [z — H^)~ x is the resolvent operator 
of H( l \ This induce, in the spirit of the Faddeev pro- 
cedure, the splitting of the wave function \^f) into two 
components 



I*) = + 1^2), 
where the components are defined by 

\i/> a )=G®vP\V), 



(10) 



(11) 



with a = 1,2. The components satisfy the set of two- 
component Faddeev-Merkuriev integral equations 



iv*> 



(12a) 
(12b) 



where Ga' is the resolvent operator of the channel 
Coulomb Hamiltonian 



H (i) = H (l) 



and the inhomogeneous term |$j ") is an eigenstate of 
H{ 1 \ 

Before going further let us examine the spectral prop- 
erties of the Hamiltonian 



(13) 



(14) 



It obviously supports infinitely many two-body chan- 
nels associated with the bound states of the attrac- 
tive Coulomb potential . The potential v$ is repul- 
sive and does not have bound states. The three-body 
potential is attractive and constructed such that 

v^\x2, 2/2) — > if ?/2 ~ * 00. Therefore, there are no two- 
body channels associated with fragmentations 2 and 3, 
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the Hamiltonian has only 1-type two-body asymp- 
totic channels. Consequently, the corresponding 



Green's operator, acting on the w^j^) term in (12a), 
will generate only 1-type two-body asymptotic channels 
in Similar analysis is valid also for \ip2}- Thus, the 

Faddeev-Merkuriev procedure results in the separation 
of the three-body wave function into components such a 
way that each component has only one type of two-body 
asymptotic channels. This is the main advantage of the 
Faddeev equations and, as this analysis shows, this prop- 
erty remains true also for attractive Coulomb potentials 
if the Merkuriev splitting is adopted. 

In the e~e~p system the particles 1 and 2, the two 
electrons, are identical and indistinguishable. Therefore, 
the Faddeev components l^i) and \ip2)-, in their own nat- 
ural Jacobi coordinates, should have the same functional 
forms 



(xiyil^i) = (2^21^2) = (xy\ip). 



(15) 



On the other hand, by interchanging the two electrons 
we have 



(16) 



where the operator V describes the exchange of particles 
1 and 2, and p = ±1 is the eigenvalue of V. Building 
this information into the formalism results the integral 
equation 



(17) 



which is alone sufficient to determine We notice that 
so far no approximation has been made, and although 
this Faddeev-Merkuriev integral equation has only one 
component, yet it gives a full account on the asymptotic 
and symmetry properties of the system. 



II. COULOMB-STURMIAN SEPARABLE 
EXPANSION APPROACH 

We solve this integral equation by applying the 
Coulomb-Sturmian separable expansion approach. This 
approach has been established in a series of papers for 
two- H and three-body [|[ [|, [| problems with Coulomb- 
like potentials. The Coulomb-Sturmian (CS) functions 
are defined by 



(r\nl) 



(n + 2Z + l)! 



1/2 



(26r)' +i exp(- 



-6r)L»+ 1 (26r), 
(18) 



with n and I being the radial and orbital angular mo- 
mentum quantum numbers, respectively, and b is the 
size parameter of the basis. The CS functions 
form a biorthonormal discrete basis in the radial two- 
body Hilbert space; the biorthogonal partner defined by 
(r\nl) = (r\nl) jr. 



Since the three-body Hilbert space is a direct product 
of two-body Hilbert spaces an appropriate basis is the 
bipolar basis, which can be defined as the angular mo- 
mentum coupled direct product of the two-body bases, 

\nvlX) a = \nl) a <g> \vX) a , (n, v = 0,1,2,...), (19) 

where \nl) a and \ v\) a are associated with the coordinates 
x a and y a , respectively. With this basis the complete- 
ness relation takes the form (with angular momentum 
summation implicitly included) 



N 

lim > 

n,y— 



nvlX) a a (nvlX\ 



lim 1 



N 



(20) 



where (xy\ni>lX) — (xy\nvlX) / {xy) . 

We make the following approximation on the integral 
equation ( p"7| ) 



\iP) = \<S>[ l) ) + G\ l >l?v\ s >pVl«\i;) 



WlJV„,(s)„ 



(s) 

i.e. the operator v\ pV is approximated in the three- 
body Hilbert space by a separable form, viz. 

1 A r 



(21) 



vfpV 



lim lfV (s Vlf 

TV^oo 



^2 MA)i v\ s > i(nVZ'A'|, (22) 



n.u.n'u'— 



where = \{nvlX\v^ pV\n' v' V X')\. Utilizing the 

properties of the exchange operator V these matrix 

Is) 

elements can be written in the form v\ = p x 

(— ) l i(nvlX\v[ s ^ \n'v'l'X')2, and can be evaluated numer- 
ically by using the transformation of the Jacobi coordi- 
nates |j. The completeness of the CS basis guarantees 
the convergence of the method with increasing N and 
angular momentum channels. 

Now, by applying the bra {n"v"l"X"\ on Eq. @ 
from left, the solution of the inhomogeneous Faddeev- 
Merkuriev equation turns into the solution of a matrix 
equation for the component vector i\) = \ (nvlX\ip) 



where 



and 



= i<ni4A|$£ ) 



(0 



nvlXlG^ln'v'l'X'}! 



The formal solution of Eq. ( p3| ) is given by 

^ = [(G«)- 1 -^ ) ]- 1 (G«)- 1 t«. 



(23) 
(24) 
(25) 
(26) 
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Unfortunately neither G\ nor $^ are known. They 
are related to the Hamiltonian H[' , which is still a com- 
plicated three-body Coulomb Hamiltonian. The approx- 
imation scheme for G^ and $j is presented in Ref. ||. 
Starting from the resolvent relation 

Gf = G 1 +G l UiG i l\ (27) 
where G\ is the resolvent operator of the Hamiltonian 

H 1= H° + uf (28) 
and the potential Ui is defined by 

U 1 =v$ ) +v£, (29) 
for the CS matrix elements (G^) -1 we get 

(G®)~ x = (Si) -1 ~U-i, (30) 

where 

G_ x = i(?i^A|G\|nWA')i (31) 

and 

U_ x = l (nvl\\U l \n'v'l'\') l . (32) 

These latter matrix elements can again be evaluated nu- 
merically. 

Similarly, also the wave function ), a scatter- 
ing eigenstate of H± , satisfies the Lippmann-Schwinger 
equation 

\$V ) ) = \$ 1 )+G 1 U 1 \S>« ) ), (33) 

where is an eigenstate of H%. The solution is given 
by 

^ } = [(G 1 )- 1 -U 1 }- 1 (G 1 )-%, (34) 

where = 1 {ni/l\\$i). 

The three-particle free Hamiltonian can be written as 
a sum of two-particle free Hamiltonians 

H° = h° Xl +h° yi . (35) 

Consequently the Hamiltonian Hi of Eq. (p8| ) appears as 
a sum of two two-body Hamiltonians acting on different 
coordinates 

Hi = h Xl + h yi , (36) 

with h Xl = h xi +Vi(xi) and h yi — h®, which, of course, 

commute. Therefore the eigenstates of Hi, in CS repre- 
sentation, are given by 

i(r^fx\^i) = i(nl\(f>i)xifx\xi), (37) 



where \<fii) and |xi) are bound and scattering eigenstates 
of h Xl and h yi , respectively. The CS matrix elements 
of the two-body bound and scattering states (nl\<fr) and 
{i>X\x), respectively, are know analytically from the two- 
body case g. 

The most crucial point in this procedure is the calcu- 
lation of the matrix elements G 1 . The Green's operator 
Gi is a resolvent of the sum of two commuting Hamilto- 
nians. Thus, according to the convolution theorem, the 
three-body Green's operator Gi equates to a convolution 
integral of two-body Green's operators, i.e. 

Gi{z) = ±-.l dz'g Xl (z-z')g yi (z'), (38) 
2m J c 

where g Xl (z) = (z - /i Xl ) _1 and g yi (z) = (z - /i Vl ) _1 . 
The contour C should be taken counterclockwise around 
the singularities of g yi such a way that g Xl is analytic on 
the domain encircled by C. 

In the time-independent scattering theory the Green's 
operator has a branch-cut singularity at scattering en- 
ergies. In our formalism Gi(E) should be understand 
as Gi(E) = lim e ^ Gi(E + is), with e > 0, and E < 0, 
since in this work we are considering scattering below the 
three-body breakup threshold. To examine the analytic 
structure of the integrand in Eq. ( ^8| ) let us take e finite. 
By doing so, the singularities of g Xl and g yi become well 
separated. In fact, g yi is a free Green's operator with 
branch-cut on the [0, oo) interval, while g Xl (E + ie — z') is 
a Coulomb Green's operator, which, as function of z\ has 
a branch-cut on the (— oo, E + ie] interval and infinitely 
many poles accumulated at E + is. Now, the branch cut 
of g yi can easily be encircled such that the singularities 
of g Xl lie outside the encircled domain (Fig. [I]). However, 
this would not be the case in the e — > limit. Therefore 
the contour C is deformed analytically such that the up- 
per part descends into the unphysical Riemann sheet of 
g yi , while the lower part of C is detoured away from the 
cut (Fig. ^). The contour in Fig. is achieved by deform- 
ing analytically the one in Fig. fi], but now, even in the 
e — > limit, the contour in Fig. avoids the singularities 
of g Xl . Thus, with the contour in Fig. the mathemat- 
ical conditions for the contour integral representation of 
Gi in Eq. (|3^) is met also for scattering-state energies. 
The matrix elements G_ 1 can be cast into the form 

Gi(E) = ^ / dz'a(E-z')g (*'), (39) 

Z7T1 J (j — 11 —V 1 

where the corresponding CS matrix elements of the two- 
body Green's operators in the integrand are known ana- 
lytically for all complex energies [0, || . 

In the three-potential formalism || the S matrix 
can be decomposed into three terms. The first one de- 
scribes a single channel Coulomb scattering, the second 
one is a multichannel two-body-type scattering due to 
the potential U, and the third one is a genuine three- 
body scattering. In our + H case the target is neutral 
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and the first term is absent. For the on-shell T matrix 
we have 



Tfi 



kfki 



w*i? (+) > + <*iS ( - 



2/ 



(40) 



where i and / refer to the initial and the final states, 
respectively, \x is the channel reduced mass and k is the 
channel wave number. Having the solutions ip an< l 

(s) 

and the matrix elements U_ l and v\ , the T matrix ele- 
ments can easily be evaluated. The spin-weighted cross 
section of the transition i — ► / is given by 



(Tfi 



na 2 (2S 12 + l)(2L+l) 2 



(2i< 



1) 



(41) 



where do is the Bohr radius, L is the total angular mo- 
mentum, 5*i2 is the total spin of the two electrons and U 
is the angular momentum of the target hydrogen atom. 



III. RESULTS 

In the numerical calculations we use atomic units (the 
mass of the electrons mi = mi = 1 and the mass of 
the proton = 1836.151527). In this paper we are 
concerned with total angular momenta L = and L = 1. 
The formula ( [io| ) gives some hint for the choice of the 
parameters in the splitting function £. We can expect 
good convergence if the "size" of corresponds to the 

"size" of "I>^ K Therefore we may need to adjust the 
parameters of the splitting function if we consider more 
and more open channels. Consequently, we also need to 
adjust the b parameter of the CS basis. We found that 
the final results and the rate of the convergence does not 
depend on the choice of b, within a rather broad interval 
around the optimal value. 

Having the T matrix we can also calculate the K ma- 
trix, whose symmetry, which is equivalent to the unitarity 
of the 5* matrix, provides a delicate and independent test 
of the method. We observed that if either the parameters 
of the splitting function are too far from the optimum or 
the convergence with the basis is not achieved the K ma- 
trix fails to be symmetric. In the separable expansion we 
take up to 9 bipolar angular momentum channels with 
CS functions up to N = 36. This requires solution of 
complex general matrix equations with maximal size of 
12321 x 12321, a problem which can even be handled on 
a workstation. We need relatively small basis because in 
this approach we approximate only short-range type po- 
tentials and the correct asymptotic is guaranteed by the 
Green's operators. 

We present first our S- wave results for energies below 
the H(n = 2) threshold. In this energy region we use 
parameters v = 2.1, xq = 3, yo = 20 and b — 0.6. Table 
| shows elastic phase shifts at several values of electron 



momenta k\. Our results, which was achieved by using 
finite proton mass, agree very well with variational cal- 
culations of Ref. 0, i?-matrix calculations of Ref. ||, 
finite-element method of Ref. [[To), as well as with the 
results of direct numerical solution of the Schrodinger 
equation of Ref. 11 1 , where infinite mass for proton were 
adopted. We also compare our calculation with the dif- 
ferential equation solution of the modified Faddeev equa- 
tions H . We can observe perfect agreements with all the 
previous calculations. 

In Table || we present S'-wave partial cross sections and 
K matrices between the H(n = 2) — H (n = 3) thresholds 
at channel energy E\ — 0.81Ry and for L — 0, where we 
have 3 open channels. We used parameters v = 2.1, xo = 
3.5, yo = 20 and b = 0.3. For comparison we also show 
the results of a configuration-space Faddeev calculation 
|l2| . We can report perfect agreements not only for the 
cross sections but also for the K matrix (except for an 
unphysical phase factor). Our cross sections are also in 
a good agreements with the results of Ref. |llj . 



In Tables II] we show the S-wave K matrix between 
the H(n = 3) — H(n = 4) thresholds at channel energy 
Ei = 0.93Ry, where we have 6 open channels. We used 
parameters v — 2.1, xq — 4, yo — 20 and b = 0.2. We 
can see that the K matrix is nearly perfectly symmetric. 
In Tables IV we present S'-wave partial cross sections be- 
tween the H(n = 3) — H(n = 4) thresholds at channel 
energies E Y = 0.93Ry, E x = 0.91Ry and E x = 0.89Ry, 
respectively. In Tables |v[- |VIH we present the correspond- 
ing P-wave K matrices and cross sections. 



IV. SUMMARY 

In this work we have studied electron-hydrogen scat- 
tering problem by solving the Faddeev-Merkuriev inte- 
gral equations. In this particular case, where two parti- 
cles are identical, the Faddeev scheme results in an one- 
component equation, which, however, gives full account 
on the asymptotic and symmetry properties of the sys- 
tem. We solved the integral equations by applying the 
Coulomb-Sturmian separable expansion method. 

We calculated S- and P-wave scattering and reaction 
cross sections for energies up to the H{n = 4) thresh- 
old. Our nearly perfectly symmetric K matrices indicate 
that in our approach all the fine details of the scattering 
processes are properly taken into account. 
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g Xi (E+iZ-z') 
— '—I 



FIG. 1: Analytic structure of g xl (E + ie — z') g yi (z') as a 
function of z' , e > 0. The Green's operator g yi (z') has a 
branch-cut on the [0, oo) interval, while g xi (E + ie — z') has a 
branch-cut on the (— oo,E + ie] interval and infinitely many 
poles accumulated at E + ie (denoted by dots). The contour 
C encircles the branch-cut of g yi . In the e — * limit the 
singularities of g xi (E + ie — z') would penetrate into the area 
covered by C . 
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g x (E+iE-z') 




FIG. 2: The contour of Fig. y is deformed analytically such 
that a part of it goes on the unphysical Riemann-sheet of 
g yi (drawn by broken line) and the other part detoured away 
from the cut. Now, the contour avoids the singularities of 
g xi (E + ie — z') even in the e — + limit. 



TABLE I: Singlet ( l S e ,p = +1) and triplet { 3 S e ,p = -1) 
phase shifts of elastic S-wave e~ + H scattering . 
k Ref. § Ref. § Ref. fTpf Ref. p| Ref. g This work 



0.1 


2.553 


2.550 


2.553 


2.555 


2.553 


2.552 


0.2 


2.0673 


2.062 


2.066 


2.066 


2.065 


2.064 


0.3 


1.6964 


1.691 


1.695 


1.695 


1.694 


1.693 


0.4 


1.4146 


1.410 


1.414 


1.415 


1.415 


1.412 


0.5 


1.202 


1.196 


1.202 


1.200 


1.200 


1.197 


0.6 


1.041 


1.035 


1.040 


1.041 


1.040 


1.037 


0.7 


0.930 


0.925 


0.930 


0.930 


0.930 


0.927 


0.8 


0.886 




0.887 


0.887 


0.885 


0.884 






3 


S e , p=-l 








0.1 


2.9388 


2.939 


2.938 


2.939 


2.939 


2.938 


0.2 


2.7171 


2.717 


2.717 


2.717 


2.717 


2.717 


0.3 


2.4996 


2.500 


2.500 


2.500 


2.499 


2.499 


0.4 


2.2938 


2.294 


2.294 


2.294 


2.294 


2.294 


0.5 


2.1046 


2.105 


2.104 


2.104 


2.105 


2.104 


0.6 


1.9329 


1.933 


1.933 


1.933 


1.933 


1.932 


0.7 


1.7797 


1.780 


1.780 


1.780 


1.779 


1.779 


0.8 


1.643 




1.645 


1.644 


1.641 


1.643 
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TABLE II: L = partial cross sections (in 7ra§) in the H(n = 
2) — H(n — 3) gap at channel energy Ei = 0.81Ry. Channel 
numbers 1, 2 and 3 refer to the channels e~ + H(ls), e~ + 
H(2s) and e~ +H(2p), respectively. For comparison the result 
of a configuration-space Faddeev calculation is presented. 

Ch.# 1 2 3 



K, 



Ki. 



1 S e , 



P 



+1 



This work 



0.564 
0.817 
0.107 



0.061 
8.373 
0.863 



0.024 
2.588 
1.722 



1.895 
-2.043 
1.798 



-2.036 
5.230 
-4.114 



1.792 
-4.114 
2.366 



Method of Ref. ]l| 



0.568 
0.814 
0.105 



0.061 
8.720 
0.824 



0.024 
2.471 
1.697 



1.864 
1.980 
-1.679 



1.971 
5.131 
-3.843 



-1.671 
-3.843 
2.028 



i S e ,p= -1 



This work 





1 


3.694 


0.001 


0.0006 


an 


2 


0.016 


10.04 


1.641 




3 


0.003 


0.547 


11.85 




1 


21.34 


0.3255 


0.6386 




2 


0.3268 


-0.4404 


-0.4161 




3 


0.6409 


-0.4161 


1.755 


Method of Ref. Jl2] 




1 


3.696 


0.001 


0.0006 




2 


0.016 


10.20 


1.678 




3 


0.003 


0.560 


11.77 




1 


24.76 


-0.3823 


-0.7510 




2 


-0.3803 


-0.4441 


-0.4167 




3 


-0.7453 


-0.4165 


1.737 
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TABLE III: S-wave K matrices in the H(n = 3) - H(n = 4) 
gap at channel energy E\ = 0.93Ry. The channel numbers 
1, 2, 3, 4, 5 and 6 refer to the channels e _ (A = 0) + H(ls), 
e"(A = 0) + H{2s), e"(A = 0) + ff(3s), e"(A = 1) + H(lp), 
e~(A = 1) + H(2p) and e"(A = 2) + fl'(ld), respectively. 



Ch.# 1 


2 3 


4 


5 


6 




Ei 


= 0.93Ry, 1 S e ,p = 


+1 




1 


1.076 


-0.647 -0.160 


0.229 


0.180 


0.074 


2 


-0.652 


1.541 -0.028 


0.129 


0.531 


0.265 


3 


-0.160 


-0.029 0.766 


0.314 


-0.757 


-0.385 


4 


0.230 


0.130 0.314 


-0.566 


-0.525 


-0.284 


5 


0.180 


0.534 -0.757 


-0.526 


0.237 


0.760 


6 


0.074 


0.266 -0.385 


-0.285 


0.760 


1.342 




Ei 


= 0.93Ry, 3 S e ,p = 


-1 




1 


9.054 


0.507 0.019 


0.666 


0.099 


0.028 


2 


0.543 


-1.700 -0.111 


-1.530 


-0.113 


-0.120 


3 


0.025 


-0.112 0.155 


-0.050 


-0.926 


-0.070 


4 


0.702 


-1.532 -0.050 


-0.851 


-0.253 


-0.048 


5 


0.104 


-0.114 -0.926 


-0.253 


0.927 


0.449 


6 


0.030 


-0.120 -0.070 


-0.049 


0.449 


-0.111 
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TABLE IV: L — partial cross sections (in nao) in the H(n = 
3) — H(n = 4) gap at channel energies Ex = 0.93Ry, Ex = 

0. 91Ry and Ex — 0.89Ry, respectively. The channel numbers 

1, 2, 3, 4, 5 and 6 refer to the channels e~(A = 0) + H(ls), 
e"(A = 0) + H(2s), e"(A = 0) + H(3s), e (X = 1) + #(lp), 
e~(A = 1) + H(2p) and e~(A = 2) + H(ld), respectively. 

Ch.# 1 2 3 4 5 6 

E x = 0.93Ry, 1 S e ,p = +1 

1 0.44 0.48(-l) 0.67(-2) 0.28(-l) 0.86(-2) 0.20(-2) 

2 0.25 3.02 0.19(-1) 0.10 0.12 0.40(-l) 

3 0.15 0.83(-l) 4.68 0.71 2.41 0.86 

4 0.49(-l) 0.34(-l) 0.55(-l) 0.49 0.59(-l) 0.24(-l) 

5 0.65(-l) 0.18 0.80 0.26 1.48 0.44 

6 0.89(-2) 0.35(-l) 0.17 0.61(-1) 0.27 2.0 
E x = 0.93Ry, 3 S e ,p= -1 

1 3.18 0.22(-2) 0.43(-4) 0.21(-2) 0.26(-4) 0.14(-5) 

2 0.12(-1) 5.92 0.93(-2) 3.77 0.44(-l) 0.61(-1) 

3 0.97(-3) 0.40(-l) 7.56 0.35 11.6 3.34 

4 0.39(-2) 1.26 0.26(-l) 0.87 O.ll(-l) 0.19(-2) 

5 0.23(-3) 0.63(-l) 3.87 0.48(-l) 9.14 1.07 

6 0.79(-5) 0.53(-l) 0.67 0.49(-2) 0.64 0.34 
Ex = 0.91Ry, 1 S e ,p = +1 

1 0.46 0.45(-l) 0.90(-2) 0.24(-l) 0.89(-2) 0.18(-2) 

2 0.26 3.74 0.24 0.77(-l) 0.74(-l) 0.15(-1) 

3 0.38 1.77 5.46 1.11 0.90 1.14 

4 0.46(-l) 0.26(-l) 0.50(-l) 0.49 0.86(-l) 0.30(-l) 

5 0.13 0.18 0.30 0.64 11.5 0.37 

6 0.16(-1) 0.23(-l) 0.23 0.13 0.22 1.15 
Ex = 0.91Ry, 3 S e ,p= -1 

1 3.26 0.22(-2) 0.23(-4) 0.20(-2) 0.20(-4) 0.17(-5) 

2 0.13(-1) 7.22 0.24(-l) 4.07 0.28(-l) 0.25(-l) 

3 0.96(-3) 0.18 9.11 0.33 0.27 12.25 

4 0.37(-2) 1.36 0.15(-1) 0.98 0.10(-1) 0.31(-2) 

5 0.26(-3) 0.69(-l) 0.89(-l) 0.74(-l) 44.97 0.44 

6 0.14(-4) 0.38(-l) 2.45 0.14(-1) 0.27 2.81 



Ex = 0.89Ry, 1 S e ,p = +1 



1 0.48 0.47(-l) 0.53(-2) 0.21(-1) 0.79(-2) 0.26(-2) 

2 0.30 4.67 0.32(-l) 0.61(-1) 0.14 0.12 

3 2.98 2.80 259.8 19.02 1.58 7.12 

4 0.45(-l) 0.20(-l) 0.72(-l) 0.51 0.77(-l) 0.13(-1) 

5 1.48 3.40 0.53 6.78 119.8 2.67 

6 0.29 2.04 1.42 0.70 1.60 38.90 



Ex = 0.91Ry, 3 S e ,p-- 



1 3.34 0.22(-2) 0.67(-5) 0.17(-2) 0.90(-5) 0.25(-5) 

2 0.13(-1) 8.68 0.94(-2) 4.33 0.17(-1) 0.87(-l) 

3 0.37(-3) 0.83 1321.0 1.75 152.8 58.26 

4 0.33(-2) 1.44 0.66(-2) 1.25 0.67(-2) 0.12(-2) 

5 0.16(-2) 0.49 50.93 0.59 124.6 56.47 

6 0.25(-3) 0.15 11.65 0.63(-l) 33.88 218.4 
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TABLE V: P-wave K matrices in the H(n = 3) — H(n — 4) gap at channel energy Ex = 0.93Ry. The channel numbers 1, 2, 
3, 4, 5, 6, 7, 8 and 9 refer to the channels e"(A = 1) + H(ls), e"(A = 1) + H(2s), e"(A = 1) + H(3s), e"(A = 0) + H(2p), 
e"(A = 0) + #(3p), e"(A = 2) + #(2p), e"(A = 2) + #(3p), e"(A = 1) + H(3d) and e"(A = 3)+H(M), respectively. 



Ch.# 1 


2 


3 


4 


5 


D 


7 


8 









= 0.93Ry, 1 S e ,p = 












1 


-1.888 


-7.518 


13.24 


9.699 


8.320 


7.148 


1.992 


4.684 


30.96 


2 


-7.525 


-29.70 


51.80 


38.14 


61. ( 6 


OQ QQ 
ZO.OO 


( .003 


1 Q OQ 


izl.o 


3 


13.30 


51.99 


-89.98 


-67.93 


-56.77 


-50.01 


-13.90 


-29.57 


-216.6 


4 


9.665 


37.98 


-67.40 


-48.21 


-42.35 


-36.39 


-9.947 


-24.07 


-156.7 


5 


8.346 


32.81 


-56.70 


-42.64 


-36.30 


-31.16 


-9.349 


-19.40 


-136.0 


6 


7.151 


28.87 


-49.82 


-36.54 


-31.08 


-28.11 


-7.718 


-17.41 


-117.3 


7 


2.006 


7.885 


-13.94 


-10.05 


-9.381 


-7.765 


-2.651 


-4.874 


-34.08 


8 


4.755 


18.64 


-29.92 


-24.51 


-19.64 


-17.67 


-4.915 


-8.953 


-73.21 


9 


31.01 


121.6 


-215.9 


-157.5 


-135.7 


-117.3 


-33.89 


-72.15 


-510.6 



Ex = 0.93Ry, 3 S e ,p = -1 



1 0.454 -0.303 -0.051 -0.020 0.080 0.043 -0.017 0.149 0.128 

2 -0.301 -2.453 -0.669 0.383 0.552 1.112 0.017 1.145 1.060 

3 -0.051 -0.672 0.398 -0.465 1.140 -0.371 0.0001 0.578 0.486 

4 -0.020 0.382 -0.464 0.354 -1.133 -0.236 0.883 -0.528 -0.110 

5 0.079 0.553 1.137 -1.136 3.936 -0.699 -3.202 1.075 -0.989 

6 0.041 1.113 -0.372 -0.236 -0.701 0.289 0.520 -0.769 -0.456 

7 -0.016 0.018 0.002 0.884 -3.203 0.518 1.673 -1.484 -0.226 

8 0.148 1.147 0.576 -0.530 1.075 -0.769 -1.483 -0.055 -0.278 

9 0.127 1.062 0.486 -0.111 -0.988 -0.457 -0.226 -0.277 0.090 



TABLE VI: L = 1 partial cross sections (in ncio unit) in the H(n = 3) — H(n — 4) gap at channel energy Ex — 0.93Ry. The 
channel numbers 1, 2, 3, 4, 5, 6, 7, 8 and 9 refer to the channels e~(A = 1) + H(ls), e~(A = 1) + H(2s), e~(A = 1) + H(3s), 
e"(A = 0) +H(2p), e"(A = 0) + H(3p), e"(A = 2) + H(2p), e"(A = 2) + H(3p), e"(A = 1) + H(3d) and e"(A = 3) + H(3d), 
respectively. 

Ch.# 1 2 3 4 5 (3 7 8 9^ 

Ex = 0.93Ry, 1 S e ,p = +1 



1 


0.380(-2) 0.104(-1) 


0.138(-2) 0.394(-l) 0.677(-2) 0.125(-1) 0.543(-2) 0.664(-2) 0.180(-2) 


2 


0.530(-l) 0.208(1) 


0.760(-2) 0.152(1) 


0.139 


0.135(1) 


0.284(-l) 


0.450 


0.103 


3 


0.319(-1) 0.321(-1) 


0.311(2) 


0.117(1) 


0.340(1) 


0.170 


0.459(1) 


0.191(1) 


0.282(1) 


4 


0.679(-l) 0.506 


0.903(-l) 0.157(1) 


0.104 


0.151 


0.213 


0.169 


0.796(-l) 


5 


0.508(-l) 0.201 


0.113(1) 


0.450 


0.415(1) 


0.103(1) 


0.169(1) 


0.113(1) 


0.871(-1) 


6 


0.219(-1) 0.448 


0.131(-1) 


i 0.150 


0.235 


0.164(1) 


0.647(-l) 0.399(-l) 0.183(-2) 


7 


0.412(-1) 0.415(-1) 


0.153(1) 


0.928 


0.169(1) 


0.282 


0.335(1) 


0.105 


0.233 


8 


0.296(-l) 0.391 


0.383 


0.440 


0.679 


0.105 


0.620(-l) 


0.800(1) 


0.283 


9 


0.807(-2) 0.890(-l) 


0.562 


0.208 


0.523(-l) 0.468(-2) 0.139 


0.283 


0.116(2) 




Ex = 


0.93Ry, 3 














1 


0.178(1) 0.484(-l) 0.853(-2) 0.158(-1) 0.603(-2) 0.167(-1) 0.457(-2) 0.191(-2) 0.625(-3) 


2 


0.247 0.235(2) 


0.390 


0.109(1) 


0.435 


0.294(1) 


0.163(1) 


0.171(1) 


0.182(1) 


3 


0.193 0.170(1) 


0.514(2) 


0.892(1) 


0.208(1) 


0.105(2) 


0.167(2) 


0.596 


0.381(1) 


4 


0.277(-l) 0.362 


0.683 


0.846 


0.576 


0.801 


0.344(-l) 0.924(-2) 0.920(-l) 


5 


0.453(-l) 0.633 


0.695 


0.251(1) 


0.373(2) 


0.525 


0.348(1) 


0.295(1) 


0.367(1) 


6 


0.291(-1) 0.981 


0.810 


0.804 


0.121 


0.416(1) 


0.887(-l) 


0.121 


0.803(-l) 


7 


0.333(-l) 0.236(1) 


0.556(1) 


0.151 


0.348(1) 


0.388 


0.276(2) 


0.290(1) 


0.186(1) 


8 


0.831(-2) 0.149(1) 


0.119 


0.240(-l) 0.177(1) 


0.315 


0.174(1) 


0.448(1) 


0.280(1) 


9 


0.259(-2) 0.158(1) 


0.760 


0.241 


0.220(1) 


0.208 


0.111(1) 


0.280(1) 


0.935(1) 
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TABLE VII: The same as in Table at channel energy Ex = 0.91Ry. 



Ch.# 1 


2 


3 


4 


5 


6 


7 


8 


9 






E x = 


0.91Ry, 1 


S e ,P = +l 












1 


0.365(-2) 


0.102(-1) 


0.102(-2) 


0.437(-l) 


0.474(-2) 


0.137(-1) 0.397(-2) 


0.521(-2) 


0.154(-2) 


2 


0.576(-l) 


0.218(1) 


0.648(-2) 


0.185(1) 


0.395 


0.122(1) 


0.694(-l) 


0.109 


0.276(-l) 


3 


0.428(-l) 


0.485(-l) 


0.963(2) 


0.213(1) 


0.380(1) 


0.262 


0.107(1) 


0.683(1) 


0.517(1) 


4 


0.829(-l) 


0.617 


0.954(-l) 


0.193(1) 


0.155 


0.245 


0.206 


0.989(-l) 


0.279(-l) 


5 


r\ c c c ( i \ 

O.ddd(-I) 


0.981 


0.127(1) 


0.115(1) 


0.509(1) 


0.605 


0.326(1) 


0.168(1) 


0.944 


6 


0.261(-1) 


0.408 


0.118(-1) 


0.246 


0.818(-1) 


0.188(1) 


0.367(-l) 


0.139 


0.520(-l) 


7 


0.563(-l) 


0.174 


0.356 


0.153(1) 


0.326(1) 


0.271 


0.570(1) 


0.157(1) 


0.262(1) 


8 


f\ A A A { 1 \ 

0.444(-l) 


0.163 


0.137(1) 


0.443 


0.101(1) 


0.616 


0.940 


0.160(2) 


0.438 


9 


0.133(-l) 


0.411(-1) 


0.103(1) 


0.125 


0.567 


0.232 


0.157(1) 


0.438 


0.187(2) 






Ex = 


0.91Ry, 3 


s e , P = -i 












1 


0.182(1) 


0.438(-l) 0.788(-2) 0.1567(-1) 0.605(-2) 


0.159(-1) 0.450(-2) 0.209(-2) 0.633(-3) 


2 


0.243 


0.241(2) 


0.209(1) 


0.170(1) 


0.126(1) 


0.479(1) 


0.106(1) 


0.632 


0.652 


3 


0.329 


0.156(2) 


0.166(3) 


0.888(1) 


0.146(1) 


0.851(1) 


0.334(1) 


0.471(1) 


0.797(1) 


4 


0.296(-l) 


0.567 


0.397 


0.139(1) 


0.464 


0.102(1) 


0.119 


0.123 


0.107 


5 


0.846(-l) 


0.311(1) 


0.486 


0.345(1) 


0.814(2) 


0.524 


0.604(1) 


0.480 


0.144(1) 


6 


0.296(-l) 


0.160(1) 


0.382 


0.102(1) 


0.706(-l) 


0.445(1) 


0.192 


0.177 


0.140 


7 


0.632(-l) 


0.263(1) 


0.111(1) 


0.885 


0.604(1) 


0.143(1) 


0.540(2) 


0.255 


0.113(1) 


8 


0.179(-1) 


0.943 


0.942 


0.550 


0.2895 


0.793 


0.153 


0.633(1) 


0.231(1) 


9 


0.551(-2) 


0.971 


0.159(1) 


0.480 


0.8620 


0.627 


0.676 


0.231(1) 


0.218(2) 



Ch.# 1 



TABLE VIII: The same as in Table at channel energy Ei = 0.89Ry. 
2 3 4 5 6 7 8 



9 



Ex = 0.89Ry, 1 S e ,p = +1 



1 0.342(-2) 0.940(-2) 

2 0.609(-l) 0.252(1) 

3 0.466 0.361(1) 

4 0.102 0.708 

5 0.495 0.205(1) 

6 0.327(-l) 0.361 

7 0.500 0.286(1) 

8 0.331 0.138(1) 

9 0.965(-l) 0.538 



0.819(-3) 0.474(-l) 0.262(-2) 0.154(-1) 0.264(-2) 0.292(-2) 0.849(-3) 

0.412(-1) 0.212(1) 0.693(-l) 0.108(1) 0.976(-l) 0.777(-l) 0.306(-l) 

0.550(3) 0.862(1) 0.384(2) 0.567(1) 0.177(3) 0.672(2) 0.232(2) 

0.326(-l) 0.233(1) 0.126 0.422 0.109 0.136 0.363(-l) 

0.128(2) 0.111(2) 0.400(3) 0.264(1) 0.166(1) 0.123(2) 0.510(2) 

0.213(-1) 0.422 0.302(-l) 0.223(1) 0.468(-l) 0.366(-l) 0.133(-1) 

0.591(2) 0.965(1) 0.166(1) 0.415(1) 0.106(3) 0.146(1) 0.291(2) 

0.135(2) 0.719(1) 0.735(1) 0.193(1) 0.879 0.695(2) 0.302(2) 

0.464(1) 0.192(1) 0.306(2) 0.706 0.175(2) 0.302(2) 0.140(3) 



Ex = 0.91Ry, s S e ,p 



-1 



1 0.186(1) 0.408(-l) 

2 0.267 0.277(2) 

3 0.376(1) 0.862(1) 

4 0.360(-l) 0.324 

5 0.108(1) 0.122(2) 

6 0.337(-l) 0.133(1) 

7 0.554 0.553(2) 

8 0.189 0.544(2) 

9 0.314(-1) 0.307(2) 



0.669(-2) 0.170(-1) 
0.983(-l) 0.972 



0.175(4) 

0.996 

0.785(2) 

0.101(1) 

0.382(2) 

0.177(2) 

0.784(2) 



0.264(3) 

0.345(1) 

0.498(2) 

0.698 

0.482(1) 

0.303(1) 

0.253(1) 



0.571(-2) 

0.415 

0.235(3) 

0.563 

0.796(3) 

0.500 

0.185(3) 

0.164(2) 

0.121(2) 



0.158(-1) 

0.398(1) 

0.266(3) 

0.698 

0.440(2) 

0.707(1) 

0.193 

0.557(1) 

0.799(1) 



0.287(-2) 

0.188(1) 

0.115(3) 

0.542(-l) 

0.186(3) 

0.227(-2) 

0.233(3) 

0.416(2) 

0.453(2) 



0.160(-2) 

0.309(1) 

0.886(2) 

0.574(-l) 

0.273(2) 

0.105 

0.694(2) 

0.106(3) 

0.981(2) 



0.258(-3) 

0.174(1) 

0.392(3) 

0.483(-l) 

0.202(2) 

0.151 

0.754(2) 

0.980(2) 

0.26(3) 



